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ABSTRACT 

We present the results of a numerical study of a galactic wind model and its 
implications on the properties of damped Lyman-a absorbers (DLAs) using cos- 
mological hydrodynamic simulations. We vary both the wind strength and the 
internal parameters of the the wind model in a series of cosmological smoothed 
particle hydrodynamics (SPH) simulations that include radiative cooling and 
heating by a UV background, star formation, and feedback from supernovae and 
galactic winds. To test our simulations, we examine the DLA 'rate-of-incidence' 
as a function of halo mass, galaxy apparent magnitude, and impact parame- 
ter. We find that the statistical distribution of DLAs does not depend on the 
exact values of internal numerical parameters that control the decoupling of hy- 
drodynamic forces when the gas is ejected from starforming regions, although 
the exact spatial distribution of neutral gas may vary for individual halos. The 
DLA rate-of-incidence in our simulations at z = 3 is dominated (80 — 90%) by 
the faint galaxies with apparent magnitude IZab < 25.5. However, interest- 
ingly in a 'strong wind' run, the differential distribution of DLA sight-lines is 
peaked at Mh a io = 10 12 h~ l M Q {IZab — 27), and the mean DLA halo mass is 
(Mdla) = 10 12 ' 4 h~ 1 M Q {IZab — 26). These mass-scales are much larger than 
those if we ignore winds, because galactic wind feedback suppresses the DLA 
cross section in low-mass halos and increases the relative contribution to the 
DLA incidence from more massive halos. The DLAs in our simulations are more 
compact than the present-day disk galaxies, and the impact parameter distribu- 
tion is very narrow unless we limit the search for the host galaxy to only bright 
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Lyman-break galaxies (LBGs). The comoving number density of DLAs is higher 
than that of LBGs down to IZab — 30 mag if the physical radius of each DLA 
is smaller than 5 h^Q kpc. We discuss conflicts between current simulations and 
observations, and potential problems with hydrodynamic simulations based on 
the cold dark matter model. 

Subject headings: quasars: absorption lines — galaxies: ISM — stars: formation 
- galaxies: evolution - galaxies: formation - methods: numerical 



Introduction 



It is widely believed that the galactic wind phenomenon has very strong effects on 
the process of galaxy f ormation. We observe winds of hot gas emanating from starburst 
galaxies such as M82 ([McCarthy et al.lll987l ). and also find evidence for la rge-scale out- 



flows in the spectrum of high -redshift starforming galaxies such as LBGs (e.g. iPettini et al. 
200ll ; lAdelberger et al.l 120031 ) . Such galactic winds heat up the intergalactic medium (IGM) 
and enrich it with heavy elements. Therefore in recent years theorists have incorporated 
phenomenological models of galactic wind into numerical simulations and studied i ts in 



fluence on galaxy formation and chemical enrichment of the IGM (ITheuns et al.l 12002 



Springel fc Hernquistll2003bl : ICen et al.ll2005l ). From those studies, it became clear that the 
observed metallicity of the Ly-a forest can be only understood if the effect of galactic winds 
is considered. In this paper, we test the model of galactic wind presently included in our 
cosmological simulations against the observations of damped Lyman-a absorbers (DLAs)Q 

Within the currently favored hierarchical A cold dark matter (CDM) model, DLAs 
observed in high redshift quasar absorption lines are considered to arise from radiatively 
cooled neutral gas in dark matter h alos. They are known to dominate the neutral hydroge n 
mass density at high redshift (e.g., lLanzetta et al.l Il995l : IStorrie-Lombardi fc Wolfel 120001 ) . 
and hence provide a significant reservoir of cold neutral gas for star formation. If this picture 
is correct, then DLAs are closely linked to star formation and must be located inside or in the 
vicinity of galaxies within dark matter halos (hereafter often just 'halos'). For these reasons, 
DLAs provide an excellent probe of high redshift galaxy formation that complements the 
study of high redshift galaxies by the direct observation of stellar light. They also provide 
useful opportunities to test hydrodynamic simulations of galaxy formation. 



J DLAs are historically defined to be absorption systems with neutral hydrogen column densities Nm > 
2 x 10 20 cm^ 2 |Wolfe et al.lll986h . 
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3), the true nature 



of DLA galaxies (i.e., galaxies that host DLAs) is still uncertain, and it is unclear how 
DLAs are distributed among dark matter halos. The obse rved large velocity widths of low 



ionization lines support a large, thick disk hypothesis (e.g. JWolfe et al.lll986l ; lTurnshek et al. 



19891 : iProchaska fc Wolfd 119971 Il998l ). while at the same time there is evidence that a wider 



range of galaxies could be DLA galaxies (e.g. 



Le Brun et al. 



1993 IK u 



Rao fc Turnsheklhoool ; Ichen fc Lanzettall2003l : IWeatherlev et alJ2005h . lHaehnelt. Steinmetz. fc Rauch 



karni et al. 



2000. 



2001 



( 119981 ) showed that their SPH simulations were able to reproduce the observed asymmetric 
profiles of low-ionization absorption lines, and argued that DLAs could be protogalactic gas 
clumps rather than well-developed massive disks that have settled down. However their 
simulations did not include the effects of energy and momentum feedback owing to star 
formation, and they only analyzed a few systems. 

To interpret the observations of DLAs, it is useful to define the 'rate-of-incidence' dis- 
tribution as a function of halo mass or associated galaxy magnitude. The 'rate-of-incidence' 
of DLAs (often just 'DLA incidence'), dN/dz, i.e., the probability of finding a DLA along 
a line-of-sight per unit redshift, could either be dominated by low- mass halos, or by very 
massive halos. In the former case, the DLA cross section in low-mass halos is significant, 
and because the halo mass function in a CDM universe is a steeply increasing function with 
decreasing halo mass, the net contribution to the DLA rate-of-incidence from low-mass ha- 
los dominates over that from massive halos. One of the critical elements that needs to be 
determined is the DLA cross section as a function of halo mass. 

Many authors have used semi-analytic models of galaxy formation based on the hierar- 
ch ical CDM model to study the distribution of DLAs. iMo et al.l (119981 ) improved the model 
of iKauffmannl (119961 ) to study the formation of disk galaxies by following the angular mo- 
mentum distribution, and examined the distribution of DLA rate-of-incidence as a function 
of halo circular velocity and impact parameter. While these calculations provide insights 
on what is expected in hierarchical CDM models, they rely on assumptions about the mass 
fraction of disks and the geometry of the gas distribution relative to those of dark matter 
halos. Furthermore they are not able to directly account for dyna mical effects from vio- 
lent merging of halos/galaxies and associated heating/cooling of gas. lHaehnelt et al.l (120001 ) 
studied the luminosity and impact parameter distribution of DLAs using simple scaling rela- 
tionships from both observations and simulations and an analytic halo mass function. There 
are some recent works that discuss the distribution and the physical properties of DLAs 



based on semi-analyt i c mo dels of galaxy formation (iMaller et al.l l200ll : lOkoshi et al.l 12004 



Qkoshi fc Nagashimal 120051 ). While these comprehensive semi-analytic models did not re- 
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quire assumptions about the cold gas mass fraction in halos, they still had to make some 
choices for the geometry of th e gas distribu t ion, s uch as an exponential radial profile of H I 
column density in the case of lOkoshi et al.l (120041 ). focusing on the virialized systems. We 
will compare our results in this paper to those of the above authors in what follows. 

On the other hand, numerical simulations are able to describe dynamical effects in a 
more realistic manner than semi-analytic models, but instead have resolution and box size 
limitations owing to finite computational resources. State-of-the-art cosmological hydro- 
dynamic simulations that evolve comoving volumes larger than ~ (10 /i -1 Mpc) 3 can now 
achieve a spatial resolution of ~ 1 kpc in comoving coordinates, but have the fundamental 
problem of being unable to produce a large population of rea listic disk galaxies with the ob- 



served number density at z = (e.g.. iRobertson et al.l 120041 ) . This is the so-called 'angular 



momentum problem' where the transfer of angular momentum from baryons to dark matter 
is excessive and the disks in the simulations become too small relative to the real ones. Even 
so, implementations of the physics of star formation and fe edback have imp roved over the 
past several years (e.g., ISpringel & Hemquistlliooi l2003al B [ Cen et al.ll2005l ). therefore it is 
of interest to compare re sults with those obtai ned from simulations that did not include star 
formation and feedback (IHaehnelt et al.lll998l ). 



In an earlier paper, iNagamine. Springel. Hernquistl (j2004bf ) studied the neutral hydro- 
gen mass density, column density distribution, DLA cross section, and the rate-of-incidence 
using a series of SPH simulations with varying box size and feedback strength. One of their 
results was that the DLA cross section in low-mass halos depends on both resolution and 
the feedback strength. This results in a strong variation in the distribution of DLA inci 



dence w hich w as often neglected in other numerical studies (IKatz et al.lll996bl ; lGardner et al. 



1997al fbL 1 20011 ). In this paper we first check that our results are not strongly affected by the 
internal parameters of the wind model that control the decoupling of the hydrodynamic 
force when the gas is ejected from a star-forming region. Then we examine the DLA rate-of- 
incidence as functions of dark matter halo mass, apparent magnitude of D LA galaxies, and 



i mpac t parameter. We compare our results with the observational results by lProchaska et al. 



( 120051 ) that are derived from the Sloan Digital Sky Survey (SDSS) Data Release 3. We have 
also improved the accuracy of the dark matter halo mass function by performing more accu- 
rate integral of the power spectrum when calculating the mass variance a(M) of the density 
field. 
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Simulations 



We use the GADGET-2 code (jSpringelll2005l ) which employs the Smoothed Particle Hy- 

drody namics (SPH) technique. It adopts the entropy-conservative formulation of lSpringel fc Hernquist 
(120021 ) which largely alleviates the overcooling problem that the previous generation of 
SPH codes suffered. Our simulati ons include radiat i ve coo ling and h eating with a uni- 
form UV backgr ound of a modified lHaardt fc Madaul (119961 ) spectrum (IKatz et al.l Il996al ; 
Dave et al.lll999l), star formation, supe rnova feedback, a phenomenological model for galac- 
tic winds (ISpringel fc Hernquistl l2003bl ) . and a sub-resolution multiphase model of the ISM 



(ISpringel fc Hernquist 



2003af ). 



We use a series of simulations of varying box size and particle number (see Table [T]) 
in order to assess the impact of numerical resolution on our results. Also, the strength of 
galactic wind feedback is varied among the 03 (no wind), P3 (weak wind), and Q3 (strong 
wind) runs, allowing us to study the consequences of feedback on our results. The adopted 
cosmological parameters of all runs are (f2 m , Q\, fib, °8, h) = (0.3, 0.7, 0.04, 0.9, 0.7). We also 
use the notation h 70 = h/0.7, where h is the Hubble parameter in units of 100 kms" 1 Mpc^ 1 . 
In this paper, we only use simulations with a box size of lOh^ 1 Mpc in order to achieve high 
spatial resolution, and focus on z = 3 since this is one of the epochs where the largest 
observational data are available and hence more accurate comparisons to observations are 
possible. 



2.1. Dependence on the wind parameters 

The details of the s tar formation model is given elsewhere ( Springel fc Hernquist 2003al lbl; 
Nagamine et al.ll2004d ). so we will not repeat them here. In the galactic wind model adopted 
in our simulations, gas particles are stochastically driven out of the dense starforming regions 
with extra momentum in random directions. The rate and amount of extra kinetic energy is 
chosen to reproduce observed mass-loads and wind-speeds. The wind mass-loss rate Mw is 
assumed to be proportional to the star formation rate, and the wind carries a fixed fraction 
X of SN energy: 



Mi 



w 



and 



(2) 



A fixed value of rj = 2 is adopted for the wind mass-loss rate, and x = 0.25 (weak wind) 
& 1.0 (strong wind) for the wind energy fraction. Solving for the wind velocity from the 
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above two equations, the two wind models correspond to speeds of v w = 242 kms -1 and 
484 kins" 1 , respectively. This wind energy is not included in esN discussed above, so in our 
simulations (1 + x) 6 sn is the total SN energy returned into the gas. 

Figure [1] compares the distribution of iVm and DLAs in a dark matter halo with a mass 
Mh ~ 1.7 x 10 12 h~ 1 M Q at z = 3 for the 03 (no wind, x = 0-0) run and Q3 (strong wind, 
X = 1.0) runs. The DLA cross section of this halo is slightly higher in the 03 run than the 
Q3 run. The DLA columns are more concentrated near the center for the 03 run than the 
Q3 run, because without the wind the neutral gas is able to cool and sink deeper into the 
potential well. In the Q3 run, the DLA columns are distributed more broadly than the 03 
run. In particular in the western side of the halo, there is a filamentary structure of DLA 
columns arising from overlapping DLA clouds. These DLA clouds may originate from the 
gas ejected by the wind but later cooled due to high gas density within the halo. 

When a gas particle goes into the wind mode, our numerical wind scheme turns off its 
hydrodynamic interactions for a brief period of time to allow the particle to escape from the 
dense star-forming region. This is done in order to obtain a well-controlled mass-flux in the 
wind, which we picture to emanate from the surface of the star-forming region. Without the 
decoupling, the mass-loading would be boosted when the accelerated wind particle kicks and 
entrains other particles in the star-forming gas. The brief decoupling is controlled by two 
internal numerical parameters that limit how long the hydrodynamic forces are ignored. The 
primary parameter Cdecoup = p/pth determines the density p in units of the threshold density 
Pth f° r star formation that the escaping wind particle needs to reach before it is allowed to 
interact normally again. The idea of letting the wind blow from the surface of the ISM is 
therefore described by the condition Cdecoup < 1. A secondary numerical parameter £ max travci 
is introduced in order to limit the maximum time t^oup = ^maxtravei/^w a wind particle may 
stay decoupled. This length scale was just added as a precaution against the unlikely - but 
in case of weak winds perhaps possible - event that a wind particle 'gets stuck' in the ISM, 
i.e. cannot climb out of the gravitational potential well far enough to reach the density 
Cdecoup Pth- We stress that our expectation is that the outcome of our simulations is rather 
insensitive to the detailed choices for these two parameters, something that we have also 
confirmed in simulations of isolated galaxies when developing the model. However, in order 
to test whether this is also the case in our much less well resolved cosmological simulations, 
we have also ran four test simulations with identical initial conditions as the original Q3 run 
but with different values of these two internal wind parameters, as summarized in Table [2j 
In particular, we have varied Cdecoup by two orders of magnitude around our default choice 
of 0.1, and we changed £ max travei between 4 and 100/i _1 kpc around our default choice of 
20/r 1 kpc. 
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The panels in Figure [2] show the Hi column density distribution and the DLA distri- 
bution for the same halo as in Fig. [I] in the four test runs. The overall distribution of gas 
is similar between the different runs, but the exact locations of DLAs in each halo are not 
identical as a result of the changes in the decoupling prescription. Some of the difference 
is presumably introduced by the randomness involved in the ejection of the gas particles 
themselves, so that statistical comparisons between the results of the different simulations 
are warranted. 

We hence examine the DLA cross section as a function of halo mass to see whether the 



distribution of DLAs is statistically different or not. iNagamine et al.l (l2004bl ) quantified the 
relationship betw een the total DLA cross section ctq° la (in units of comoving /z~ 2 kpc 2 ; note 



that Fig. 2 — 4 of INagamine et al.l ( 2004bl ) plotted comoving /i 70 2 kpc 2 ) and the dark matter 

log (T^LA = « (log M halo - 12) + )3, (3) 



halo mass fin units of h 1 M m ) at z = 3 as 



with slopes a = 0.72, 0.79, 0.84, 0.93, 1.02 and the normalization (3 = 3.94, 3.99, 3.98, 4.03, 4.18 
for the 03, P3, Q3, Q4, and Q5 runs. The slope a is always positive, and the massive halos 
have larger DLA cross section, but they are more scarce compared to less massive halos. The 
quantity (3 gives the value of logoff at M^aio = 10 12 h~ l M Q . Two qualitative trends were 
noted: (1) As the strength of galactic wind feedback increases (from 03 to Q3 run), the slope 
a becomes steeper while the normalization (3 remains roughly constant. This is because a 
stronger wind reduces the gas in low-mass halos at a higher rate by ejecting the gas out of 
the potential well of the halo. (2) As the numerical resolution is improved (from Q3 to Q5 
run), both the slope and the normalization increase. This is because with higher resolution, 
star formation in low-mass halos can be described better and as a result the neutral gas 
content is decreased due to winds. On the other hand, a lower resolution run misses the 
early generation of halos and the neutral gas in them. 

Figure [3] shows the distribution of DLA cross section as a function of dark matter halo 
mass for the four test runs with different decoupling parameters. Open triangles are the 
median in each halo mass bin. The solid line is the power-law fit to these median points (see 
Table EJ), and the dashed line is the fit to the original Q3 run. The shaded contours in the 
background give the actual distribution of halos equally spaced in logarithmic scale. Table [3] 
shows that the distributions of DLA cross sections in the four test runs are rather similar 
to that of the original Q3 run. There is a slight tendency that the slope is shallower in 
the four test runs than in the original Q3 run, something that could be related to extensive 
changes in the time integration and force calculation scheme of the simulation code relative 
to the older code version used in the original Q3 run, rather than to systematics introduced 
by the decoupling scheme. In any case, the fluctuations in our results introduced by drastic 



- 8- 



variations of the internal numerical parameters used in the decoupling scheme are much less 
than the large differences between the Q3 (strong wind) and 03 runs (no wind), meaning 
that the decoupling parameters are not a significant source of uncertainty in our results. 

Figure H] shows the Hi column d ensity distribution funct ion and the cumulative rate-of- 



incidence for the four test runs (see iNagamine et al.l l2004bl . for the definition of these two 



quantities). These figures also show that the statistical distribution of DLAs in the four test 
runs is similar to the original Q3 run. Likewise, the impact parameter distribution which 
we will discuss in Section E] is also very similar to the original Q3 run for the four test runs. 
Therefore we focus on the original 03, P3, Q3, and Q5 runs in the following sections. 



3. Differential distribution of rate-of-incidence as a function of halo mass 

The differential distribution function of DLA incidence can be computed as 



diVi 



DLA 



^[Mn(M,z) \n(lO)} a c D ° LA (M,z) 



(4) 



c/H(z) with H (z) 



dz d log M dz 

where n(M,z) is the dark matter halo mass function, and dr/dz 
H E(z) = H ^/n m (l + z) 3 + tt A for a flat universe. We use the Isheth fc Tormenl (jl999l ) 
parameterization for n(M) as shown in Figure [5j Note that the dependence on the Hubble 
constant disappears on the right-hand-side of Equation (jlj) because dr/dz scales as h -1 , 
Mn(M) scales as h 3 , and o"dla scales as h~ 2 in the simulation. For t he cum ulative version 
of this calculation, see Equation (8) and Figure 5 of INagamine et al.l (]2004bl ). Equation (Tj0) 
can be derived from the following expression for the DLA area covering fraction on the sky 
along the line element c dt: 



cL/Vdla = n phys {M)dM -a^l -cdt 



(l + zfn co (M) dM-o 
n co {M) dM-a^ LA -dr, 



phys 
DLA 



a dr 



(5) 
(6) 
(7) 



where we have used cdt = adr and ct^la = (1 + z ) "dla- Here a is the scale factor, 
dr is the line element in comoving coordinate, and n co (M)dM and n p ^ ys (M)dM are the 
comoving and physical number density of halos in the mass range [M, M + dM], respectively. 
Sometimes the 'absorption distance' dX is defined as dX = — (I + 2;) 3 cdt = — (1 + z) 2 dr = 
jjfc (1 + z) 2 dz = (1 + z) 2 dz/E(z), and is used to express the rate-of-incidence as 



(8) 
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For z = 3 and our adopted cosmology, dX/dz = 3.5867. In Equations (jSJ) to (jSJ), we left in 
the dependence on halo masses explicitly, but in practice an integral over a certain range of 
halo mass has to be performed when comparing with actual observations. 

We now use the power-law fits for crg 3 LA (M, z) described above to compute the differential 
distribution of DLA incidence via Equation (111) . The result is shown in Figure [6] for all 
the simulations at z = 3. The qualitative features of the curves are easy to understand. 
Because n(M) oc M~ 2 at M « 10 8 — 10 12 h~ l M & (see Figure E]), the distribution is flat when 
c"dla oc M. In fact, n(M) is slightly shallower than M -2 (more like M -1,8 ), therefore the 
distribution for the P3 run is almost flat at 10 s < M < 10 12 /i _1 M , because <tdla oc M a79 
in this simulation. At masses higher than 10 12 h~ 1 M Q , the mass function deviates from the 
M~ 2 power-law significantly, and the distributions for all runs quickly drop off to a small 
value. 

The halo masses where each distribution peaks are listed in the second column of Table HI 
The peak halo mass M pea k becomes larger as the feedback strength increases. For the 03 
run, we indicated M pea k = 10 8,5 h~ x M & in parentheses because we think that the DLA cross 



section rapidly falls off at this halo mass based on the work by iNagamine et al.l (l2004bl ) and 
the peak halo mass is simply this cutoff mass-scale. The peak halo mass is significantly 
larger for the Q4 (M pcak = 10 116 ^ _1 Af ) and Q5 (M peak = 10 12 Ii^Mq) runs compared to 
other runs. 



4. Mean & Median halo masses of DLAs 

For each distribution shown in Figure EJ we compute the mean DLA halo mass as 

(M dla ) - r div dl m~ (9) 

JO dzdlogA/ U1U & iW 

J oo M 2 n(M)a DLA (M)dlogM 



f™Mn(M) a DLA (M) dlogM 



and the result is summarized in TableHl The mean halo mass is smaller for the 'no- wind' (03) 
run, and is larger for the 'strong-wind' (Q3 to Q5) runs. This is because of the steepening of 
the relationship between odla and M halo as the feedback strength increases. But the mean 
halo mass is in the range Mdla — 10 11 ' 5 — 10 12 ' 5 /i _1 M Q for all runs. While this mass-scale 
is not as large as a Milky Way sized halo (Mh a i — 10 12 — 1O 13 M ), it is certainly more 
massive than th at of dwarf galaxie s. We also list the mean of log Mdla for comparison with 



the results by iBouche et al.l (120051 ) . We will discuss the implications of these mass-scales in 
Section [HI 
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One can also look at the median halo mass M50, below (or above) which 50 percent of 
the DLA incidence is covered. The quantity M 50 is defined by the following equation: 

/*oo /»oo 

0.5=/ n(M) a DLA {M) dM / I n(M) a DLA (M) dM. (11) 
Jm 50 Jo 

The median halo mass is always smaller than (Mdla) 101 all runs. The values of M 50 are 
indicated by the open crosses in Figure El and (M DL a) are indicated by the filled squares. 
We emphasize that M 50 is much smaller than (Mdla), and that (Mdla) could be biased 
towards a large value because of the weighting by the halo mass. 

Furthermore, in Table H] we also give the 75 percentile halo mass M75, below which 75 
percent of the DLA incidence is covered. For example, in the Q5 run, halos with masses 
M ha i < 10 12 ' 3 /i _1 M are responsible for 75% of the DLA incidence at z — 3. The values of 
M75 are shown with open triangles in Figure El 



5. Luminosity distribution of DLA galaxies 



The luminosity distribution of DLA gala xies constrains their nat ure, facilitating the in- 
terpretation of observations of DLA galaxies. iNagamine et al.1 (J2004eJ) computed the spectra 
of galaxies in the same sim ulations used in this paper using the population synthesis model of 
Bruzual fc Charlotl (120031 . hereafter BClib03) based on the stellar mass, formation time, and 
metallicity of each stellar particle that makes up simulated galaxies. Using the computed 
spectra, we sum up the monochromatic luminosity at rest-frame 1655A (chosen because it 
corresponds to the observed-frame 71 band of the U n G7Z system) of all the galaxies enclosed 
in the max imum radius of ea ch dark matter halo identified by the f riends-of-friends g rouping 
algorithm (IDavis et al.lll985l ). The following formula (see Eq.(2) of lNight et aDbood ) is used 
to compute the absolute AB magnitude at 1655A, because BClib03 outputs its spectrum L\ 
in units of L & A" 1 where L & = 3.826 x 10 33 ergs _1 : 



M AB = -2.51og(A 2 L A ) + 13.83. 



(12) 



Figure [7] shows the relationship between M AB and the halo mass for 03, P3, Q3, and Q5 
runs at z = 3. The three dashed lines correspond to 



M 



AB 



^(logMh^-l^-d, 



(13) 



where C\ = 23.5,22.5, and 21.5 from top to bottom, and Mh a io is in units of Ji^Mq. We 
adopt C\ = 23.5 for the 03 run, C\ = 22.5 for the P3 run, and C\ = 21.5 for the Q3 and Q5 
runs. This figure shows that, on average, the galaxies in the 'no-wind' run (03) are brighter 
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than those in the ' strong-wind' run (Q3 &: Q5) by about a factor of six. This was pointed 
out in Figure 6 of iNagamine et al.l (j2004el ) and was attributed to the suppression of star 
formation by a strong galactic wind. 



Here, the apparent and absolute magnitudes are related to each other as 
m AB = M AB + 2.5 log(l + z) + 5 log(d L /10pc), 



(14) 



where di is the luminosity distance. Notice the positive sign in front of the 2nd term (which 
is normally neg ative); this is beca use of the definition of M AB in Equation ( fl2l) (see also 
Equation (2) of iNight et al.ll2006l ). Inserting Equation (fT3|) into Equation (TTJJ and using 
d L (z = 3) = 2.542 x lO^^Mpc for our flat A cosmology, we obtain the relation between 
the apparent magnitude TZab and the halo mass Mh a io as 



TZ AB = -2.5 log M halo + C 2 - 5 log h 



70- 



(15) 



where C 2 = 55.03 (03 run), 56.03 (P3 run), and 57.03 (Q3 and Q5 run), and Mha\n is in units 
of h~ l M Q . These numerical values are consistent with the ones adopted by lHaehnelt et al. 



( 120001 ) . where they assumed m AB = 26.6 — 7.5 log(t> c /200kms x ) = —2.5 logM halo + 55.7 for 
a flat A cosmology. The circular velocity v c at a radius of overdensity 200 is computed as 
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(16) 
(17) 



for our flat A cosmology and p is the mean density of the universe at redshift z. 



Figure [8] shows the cumulative distribution for the absolute value of DLA rate-of- 
incidence as a function of ap parent H-ab magnitude o f host galaxies. The slight differences 
in these results from Fig. 5 of INagamine et al.l (l2004bl ) are due to an improved calculation of 
the halo mass function with a more accurate integration of the power spectrum when com- 
puting the mass variance a(M). The yellow hatched region sho ws the observational est imate 
log(dAYd;z) = -0.60 ± 0.10 at z = 3 obtained from Fig. 8 of IProchaska et al.l J2005h . We 
estimate that the range of dN/dX centered at z = 3 is [0.065, 0.09] with a central value 0.07, 
and then multiply dN/dX by dX/dz = 3.587 for z = 3 to obtain the above observational 
value. We find that the P3, Q3 and Q5 runs underpredict dN/dz. We will discuss this 
discrepancy further in Section [HI 

It is clear from Figure [8] that only a very small fraction of DLA incidence is associ- 
ated with galaxies brighter than TZ A b — 25.5. This is more evident in Figure [9] where the 
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cumulative probability distribution of DLA incidence is shown. In panel (a), we normal- 
ize the cumulative distribution by the value at M halo = 10 98 h~ 1 M Q , or equivalently, halos 
less massive than the above value are assum ed not to host DLAs at all. Here we roughly 
reproduce the result of lHaehnelt et al.l (120001 ) where they adopted the cutoff circular veloc- 
ity of 50kms _1 , which corresponds to Mh a io — 10 9 ' 8 h~ l M & . Similarly to their result, we 
find that for this cutoff mass, only 10 — 20% of DLA sight-lines are contributed by galaxies 
brighter than the spectroscopic limit TZab = 25.5, and 70 — 90% of the DLA sight-lines are 
contributed by galaxies brighter than magnitude TZab = 30. 



However, according to the work by iNagamine et al.l (l2004bl ). only halos with masses 



Miaio > 10 8 ' 5 h 1 M (or equivalently v c > 18kms x ) contribute to the DLA cross s ection 



(see their Figure 2 & 3), so the cutoff velocity of 50kms _1 adopted by lHaehnelt et al.l (120001 ) 
might be too high. In Figure [Hb we show a case with cutoff halo mass Mh a i > 10 8 ' 5 h~ 1 M & . 
Here, less than 15% of DLAs are contributed by galaxies with IZab < 25.5, and 50 — 70% 
by those with TZab < 30. 



6. Impact parameter distribution 



We also compute the impact parameter distribution (i.e., pr ojected separatio n betwe en 
DLA sight-lines and the nearest galaxy) using our simulations. INagamine et al. (h004bl lc}) 
calculated the neutral hydrogen column density Ahi for each pixel with a size of e 2 on pro- 
jected planes covering the face of each dark matter halo, where e is the smoothing length of 
the simulation. The pixels with Nm > 2 x 10 20 cm -2 were equally regarded as DLA sight- 
lines. Therefore there are multiple DLA sig ht-lines per halo . Galax ies in our simulations 
are identified as collections of star particles. Nagamine et al. ( 2004el ld1) computed the spec- 
trophotometric properties of simulated galaxies based on the formation time, stellar mass, 
and metallicity of individual star particles, and showed that the simulated luminosity func- 
tions agree reasonably well with the observations if a mean extinction of E(B — V) = 0.15 
is assumed. This level of extinction is co nsistent with obsery ationally derived values for the 
Lyman break galaxies (LBGs) at z ~ 3 (jShapley et al.ll200ll ). 



Knowing the locations of both DLAs and galaxies in the simulation, we compute the 
impact parameter for each DLA pixel by searching for the nearest galaxy on the projected 
plane. If a nearby galaxy cannot be found within the same halo, we allow the search to 
extend further. Figure [10] shows the cumulative probability distribution of DLA incidence 
as a function of impact parameter 6 p h ys (in units of physical kpc). Figure ITUa shows the 
results of the 03, P3, and Q3 runs to highlight the impact of galactic wind feedback. As the 
feedback strength increases, gas in low mass halos is ejected more efficiently, and the neutral 
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hydrogen content decreases. Therefore, the relative contribution from higher mass halos 
increases and the impact parameter distribution becomes broader. Another notable feature 
of this plot is that, if all galaxies are allowed to be a candidate DLA galaxy no matter how 
faint they are, then the majority (over 90% for the 03 run, and 80% for the P3 and Q3 runs) 
of DLA sight-lines have the nearest galaxy within & p h ys = 5 h^Q kpc. However, if we limit the 
search for the nearest galaxy to those brighter than TZab = 30 or 28 mag, then a large fraction 
of DLAs, in particular those in low mass halos, cannot be associated with a qualified galaxy 
within the same halo, resulting in a much broader impact parameter distribution. About 
30% of all DLA sight-lines in the 03 run have impact parameters 6 p h ys > 5 kpc for the 
limited search of a nearest galaxy with TZab < 28. 

Figure [TUb shows the result of the Q3 and Q5 runs. The higher resolution run (Q5) has 
a narrower impact parameter distribution than the lower resolution run (Q3), because it can 
resolve more low-mass galaxies which host DLAs with low impact parameters. Therefore the 
relative contribution from DLAs with low impact parameters becomes larger. Also, galaxies 
in massive halos will be better resolved in the Q5 run than in the Q3 run, and this will also 
result in smaller impact parameters for the DLAs in massive halos. We expect that, in a 
run with even higher resolution, the impact parameter distribution will remain as narrow as 
that of the Q5 run, because increasing the resolution always seems to increase the relative 
number of columns with low impact parameters. We also show the case where we limit the 
search for the nearest galaxy to those brighter than TZab = 30 mag. In this case, results 
from the Q3 and Q5 runs agree well, and about 40% of all DLA sight-lines have impact 
parameters 6 p h ys > 5 kpc. 

For most of the DLAs with low impact parameters, there is a galaxy within the same 
dark matter halo. However, for a few to 10% of the DLA sight-lines that are in the low mass 
halos in 03, P3, and Q3 runs, there are no galaxies within the same halo and the nearest 
galaxy is in another halo, as indicated by the offset of the curves at the bottom right corner 
of the plot in panel (a). In the Q5 run, there are almost no DLA sight-lines that do not have 
a galaxy within the same halo, and the distribution approaches zero at large & p h ys values in 
panel (b). 

Over all, our impact param eter distribution seems to be much narrower than that ob- 



ta ined bvlHaehne 



of Haehnelt et al 



t et al.l (120001 ) if we do not limit our search to the bright galaxies. Figure 3 
( 200ol ) suggests that 60% of DLA sight-lines have b < 1 arc second, but 



our results indicate that more than 80% of DLAs have b < 1 arc second if we do not limit 
our search for the nearest galaxy to those with TZab < 30 mag. The differences between 
the two results might stem from subhalos within the massive halos in the simulations , which 



Haehnelt et al.l (120001 ) did not take into account. The calculation of lHaehnelt et al.l (120001 ) 
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assigns a single effective DLA radius to each halo and assumes that the DLAs cover a circu- 
lar area centered on each halo, whereas in our simulations, massive halos contain numerous 
galaxies and the geometry of the DLA cross section cannot be characterized by a circular 
area centered on each halo. For example, the most massive halo in the 03 run contains 
143 galaxies, and in the Q5 run 1110 galaxies. Therefore the DLAs in the simulation would 
be able to fin d the nearest ga l axy a t distances much smaller than the effective DLA radius 
computed by lHaehnelt et al.l (120001 ) . But when we limit our search for the nearest galaxy 
to those with TZa p < 30 mag, then ou r results for the Q3 and Q5 run become similar to 
that obtained by lHaehnelt et al.l (120001 ) . It is not very clear how our results compare to 
Figure lid by lMo et al.l (119981 ). as these authors examined only the differential distribution 
of impact parameter and not the cumulative probability distribution. But their differential 
distribution peaks at 2.5 h~ l kpc and a significant fraction (more than 50%) of the area under 
the curve appears to come from b < 5 h^ 1 kpc, in reasonable agreement with our results. 



7. Number density of DLAs 

Finally, we discuss the comoving number density of DLAs. Assuming that the charac- 
teristic covering area of each DLA A DL a is fixed with a physical radius r DL A (i-e., ^4 DL a = 
7r 7"j5 LA ) , we can compute the cumulative number density of DLAs as a function of halo mass 
as 



poo ff phys 

N BhA (>M)= dMn(M)-^A, (18) 

JM ^DLA 

where the definitions of n{M) and (7^1 are the same as described in Section [3j Note 
again that odla is the total DLA cross section of each dark matter halo; in other words, 
if there are 100 DLA gas clouds in a massive halo, then this halo has a total DLA cross 
section of cxJ^a — 100 ^dla- Then, using Equation (TTSj) allows us to obtain the cumulative 
comoving number density of DLAs as a function of apparent TZab magnitude A DL a(< TIab) 
as shown in Figure [Hj Here, three different values of physical radius for DLAs are assumed: 
^dla = 1,5 and 20 kpc. For each radius, the results of four simulations (03, P3, Q3, and 
Q5 runs) are shown. Also given is the cumulative number de nsity of LBGs, At ; rq ( < m) , 



computed by integrating the Schechter luminosity function of lAdelberger fc Steidell (120001 ) 
with parameters (m*, a, $*[/i 3 Mpc^ 3 ]) = (24.54,-1.57,4.4 x 10~ 3 ): 

A LBG (< m)= I $(m)dm, (19) 

J m 

where $(m) = (0.4 In 10) $* 10^ (a+1) exp(-lO^) as a function of apparent magnitude m, 
and jj, = 0.4 (m* — m). The comoving number density of LBGs at z — 3, A^bc = 4 x 10 -3 
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( lAdelberger et al.ll2003l ). is also shown as the data point at the magnitude limit of IZab = 
25.5. 

Figure [TT] shows that the comoving number density of DLAs is larger than that of 
LBGs down to the magnitude limit of IZab — 25.5 if the physical radius of each DLA is 
smaller than tdla — 5 kpc. The two number densiti e s roug hly agree with each other 
at IZab = 25.5 when tdla — 20 Ii^q kpc. Earlier, ISchayd (j200ll ) argued that the observed 
DLA rate-of-incidence can be accounted for if each LBG were accompanied by a DLA cross 
section of irr 2 = Trf!^ 1 kpc) 2 assuming (dN/dz) DLA = 0.20 and JV LB g = 0.016 h 3 Mpc~ 3 
(down to IZab = 27 mag). His result is in good agreement with the result of the Q5 run 
with ?"dla = 20 kpc shown in Figure [TTJ However, this large DLA radius is somewhat 
unrealistic because this model implies that all halos with masses Mh a i > 10 12 /i _1 M Q host 
such a large disk at z = 3 that are responsible for DLAs, and none of the less massive halos 
host DLAs at all. This picture is quite the contrary to our simulation results that indicate 
halos down to masses Mh a i = 10 8 ' 5 h~ l M Q could host DLAs, and simulated DLAs in lower 
mass halos are clumpy and smaller in cross section. 

The existence of extended disks at high redshift can be observationally tested by search- 
i ng for extended emission from stars in de ep imaging data such as the Hubble D eep Fields 
( Bouwens et al.ll2003l ; iFerguson et al.ll2004l ) . The most recent study on this issue by lBouwens et al 
(120041 ) using the Hubble Ultra Deep Field (HUDF) suggests that high redshift U BVi-dxoY>ovX 
galaxies are compact in size (~ 0.1 — 0.3 arc seconds) and that extended sources (> 0.4 arc 
sec , > 3 kpc) are rare. Another observational analysis of high redshift galaxies in HUDF 
by IWolfe fc Chen! (120061 ) also suggests that there are no extended disks down to very faint 
surface brightness. Our simulation results and impact parameter distribution are in accord 
with these observational results. If our picture is realistic and each DLA has a physical size 
of < 5 kpc, then it means that there are multiple clumps of DLAs around massive galaxies 
such as LBGs at z ~ 3, although the fraction of DLA incidence covered by the DLAs in 
massive halos (Mh a i > 10 12 h~ l M Q ) is smaller than that in less massive halos that host very 
faint galaxies {IZab > 27 mag). 



8. Discussion & Conclusions 

Using state-of-the-art cosmological SPH simulations, we performed a numerical study 
of a galactic wind model, and examined the distribution of DLA rate-of-incidence as a 
function of halo mass, galaxy apparent magnitude, and impact parameter. We find that 
the majority of DLA rate-of-incidence in the simulations is dominated by relatively lower 
mass halos (M ha i < 10 12 /i~ 1 M Q ) and faint (JZab > 25.5) galaxies. This conclusion agrees 
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with the generic prediction of the semi-analytic model of iKauffmannl (119961 ). that the DLA 
galaxies at high redshift will typically be smaller, more compact, and less luminous than disk 
galaxies at the prese nt epoch, although th i s anal ysis was restricted to an Q m = 1 universe. 
More recent work by lOkoshi fc Nagashimal (120051 ) also suggests that the low- redshift (z < 1) 
DLA galaxies are mainly low-luminosity, compact galaxies. Combined with our results, the 
dominance of faint galaxies in DLA incidence seems to be a generic prediction of a CDM 
model. This can be ascribed to the steeply rising dark matter halo mass function towards 
lower masses in CDM models, and to the fact that the small DLA cross sections in these low 
mass halos add up to a large portion of the total DLA incidence when multiplied by a large 
number of low- mass halos. 



8.1. On the DLA Halo Mass 

We characterize the differential distribution of DLA rate-of- incidence dN/ (dz d log M) 
with various halo masses listed in Table HI We find that the mean DLA halo mass increases 
with increasing galactic wind feedback strength, because winds are able to eject the gas in 
lower mass halos, suppressing their D LA cross section, result ing in a larger relative contribu- 
tion from higher mass halos (see also Nagamine et al. 2004bl ). The mean DLA halo mass for 



the Q5 run was found to be (M DLA ) = 10 12 - 4 h~ l M & and (log M DLA [fc _1 M Q ]) = 11.3 when we 
limit the DLA distribution to Mh a i > 10 8 ' 5 /i _1 M Q . The latter value is close to that obtained 



by lBouche et al.l (120051 . (log Mdla[ ^Mq]) = 10.9), but this comparison is not fully appro- 
priate because their simulation only resolved halos with masses log(Mh a io[ h^Mo}) > 10.5 
and they did not attempt to extrapolate their halo distribution using an analytic halo mass 
function. Therefore, the value of mean DLA halo mass they derived was biased toward a 
larger value owing to limited resolution, and it is an upper limit as the authors described 
in their paper. Since their simulation did not include galactic wind feedback, their result 
should be compared to that of our 03 ('no- wind') run. In fact, the value of concern for the 
03 run is (logMoLAf^ -1 ^©]) = 10.4, and it is lower than their halo mass resolution when 
we take the full distribution of DLAs into account down to M halo = 10 8 ' 5 h~ l M & We also 
note that there are large differences between (logMoLA), log(ikfrjLA), & therefore the 
former two quantities are not to be confused with M 50 . 

Given that the results of the Q3 and Q5 runs do not agree, the results given here may 
ap pear as if they haye not c onverged yet. According to the results of the R-Series presented 



by iNagamine et al.l (l2004bl ). we consider that the lowest halo mass that could host DLAs 
is Mhaio ~ 10 8 ' 5 /i _1 M Q . In the Q5 run, there are 150 dark matter particles for a halo of 
this mass, therefore it is close to the resolution limit. The Q5 run should be close to the 
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convergence if not fully converged. It is possible that in a run with higher resolution run 
than the Q5 run, the mean DLA halo mass could be even higher than that of the Q5 run. 

There have already been several observatio nal attempts to con s train DLA halo masses 



via cross-correlation between D LAs and LBGs (IGawiser et al.ll200ll ; lAdelberger et all 12003 



Bouche fc Lowenthall 120031 . 12004 ). b ut owing to l i mited sample sizes, the results have been 
mostly inconclusive. More recently, ICooke et al.l (120061 ) has measured the cross-correlation 
between 11 DLAs and 211 LBGs, and constrained the DLA halo mass to be ~ 10 11,2 M Q . It 
is encouraging tha t this measurement is close to our results listed in Table HI The result of 
Cooke et al.l (120061 ) suggests that at least some DLAs are associated with relatively massive 
halos, close to the LBG halo masses (~ 1O 12 M ). There could of course be some distribution 
in the LBG halo masses (perhaps from M ha i = 10 10 to 1O 13 M ; see the broad distributio n 
of stellar masses at the magnitude limit TZab = 25.5 in Figure 4 of iNagamine et all (j2004el )). 



and likewise DLAs could also have a broad halo mass distribution. The majority of DLA 
sight-lines in the simulations are dominated by lower mass halos in spite of a relatively large 
mean DLA halo mass, which is also reflected in the large differences between the median 
halo mass and the mean halo mass as we summarize in Table HI 



8.2. On the Luminosity Distribution 

As for the luminosity distribution of DLA galaxies, we find that only about 10% of DLA 
sight-lines are associated with galaxies brighter than TZab = 25.5 mag. This suggests that 
only about 10% of DLA galaxies will be found in searches for the bright LBGs at z ~ 3 
down to TZab — 25.5. The dominance of DLA sight-lines by the faint galaxies is a generic 
result in all of our simulations, therefore we consider that this conclusion will not change in 
a higher resolution run than the Q5 run. 



We reproduce the result of lHaehnelt et al.l (120001 ) when we cut off the DLA distribution 



at M halo = 10 9 8 /i _1 M Q (or equivalently a circular velocity of v c = 50 km s -1 ), and in this case 
70 — 90% of DLA sight-lines are associated with galaxies brighter than TZab = 30 mag. This 
agreement with the Haehnelt et al. result is not surprising, because they assumed a relation 
&dla oc v * oc with a ~ 2.5 — 3 and a normalization matched to the observed dN/dz of 

DLAs, and our SPH simulations suggest a = 2.2 — 3.1 depending on the feedback strength 
and resolution. However, our simulations indicate that even lower mass halos contribute to 
the DLA cross section, and when we cut off our DLA distribution at Mh a i = 10 8 ' 5 h~ 1 M Q , 
only 50 — 70% of DLAs are associated with galaxies brighter than TZab = 30 mag. Thus, 
it would be possible to detect DLA galaxies with ~ 50% efficiency by searching down to 
TZab = 30 mag. 



-18- 



While the dominance of faint galaxies among DLAs seems to be a generic prediction of 
the CDM model, we note that the recent photometric survey of low-redshift (z < 1) DLA 
galaxies by IChen fc Lanzettal (120031 ) on the contrary suggests that a large contribution from 
dwarf galaxies is not necessary to account for their observed DLA inc idence. However, this 
conclu sion might be affected by the 'masking effect' emphasized by lOkoshi fc Nagashima 
( 120051 ). Current DLA surveys might be biased against the detection of DLAs associated 
with faint and compact galaxies, because such galaxies would be buried under the bright 
QSO that has to be masked for the detection of a DLA galaxy. This effect would be more 
severe if the impact parameters are small as our present work suggests, as well as some 
observational studies of high redshift D LA galaxies that imply very small impact parameters 
less than se v eral kpc (i.e., < 1 a rcsec) iFvnb o et al.lll999l ; iKulkarni et al.ll2000t iMoller et al.l 



20021 . 12004 ). iFynbo et al.l (Il999l ) suggests that, based on the properties of a limited sample 
of observed high redshift candidate DLA galaxies, a large fraction (~ 70%) of DLA galaxies 
at z ~ 3 could be fainter than TZab = 28 mag. It is possible that ev olution with redshift 



i s rela tively strong as suggested by the chemical evolution model of lLanfranchi fc Friaca 



( 120031 ) . in the sense that the high redshift DLAs are dominated by dwarf galaxies and 
the low-redshi ft ones by disks. But the la tter scenario seems to be inconsistent with the 
predictions by lOkoshi &: Nagashimal (120051 ). 



Hopkins et al.l (120051 ) also argued for the dominance of faint galaxies for the DLA galax- 



ies based on the comparison of global quantities such as the density of gas mass, stellar mass, 
metal mass, and star formation rate. They also suggested that the DLAs may be a distinct 
population from LBGs, but we note that the dominance of faint galaxies for the DLA rate- 
of-incidence does not immediately mean that DLAs do not exist in LBGs. It simply means 
that the area covered by the DLAs associated with LBGs is much smaller than that in faint 
galaxies. We plan to investigate the connection between DLAs and LBGs in more detail in 
future work. 



8.3. On the Impact Parameter Distribution 



Our distribution of impact parameters is significantly narrower than that of lHaehnelt et al. 



( 120001 ). and we ascribe this difference to substruct u res w ithin the massive halos in our sim- 
ulations. The differential distribution of iMo et al.l (119981 . Figure lid) seems to be in better 
agreement with our results, with roughly half of DLAs having 6 phys < 3/i -1 kpc, but this com- 
parison is probably not fully appropriat e beca use their model is restricted to centrifugally 
supported disks. Indeed, lHaehnelt et al.l (120001 ) comment that their effective DLA radius is 
about a factor of 10 larger than the expected scale length of a centrifugally supported disk 
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if the angular momentum of the gas owes to tidal torquing during the collapse of a dark 
matter halo. However, we caution that the large effective DLA radius for a massive halo 
does not necessarily mean that DLAs consist of only large disks centered on halos; i.e., the 
DLA cross section could be distributed to hundreds of galaxies (associated with subhalos) 
embedded in a mas sive halo, with e ach D LA gas cloud being fairly compact. This point was 
also emphasized by iGardner et all (120011 . Section 4.1). 



Ou r narrow impa c t para meter distribution at first sight might seem inco nsistent with 
that of IGardner et all (120011 ). but in fact they are not. I Gardner et al.1 (120011 ) mostly used 
SPH simulations with 64 3 particles which were able to resolve halos only down to Mh a io = 
8.2 x 1O 1O M , and found that nearly all DLAs lie within 15 — 20 kpc of a galaxy center. Given 
their mass resolution, they were not able to simulate galaxies fainter than TZab — 30 mag at 
z = 3. If we limit the search for the nearest galaxy to those brighter than TZab — 30 mag (see 
Figure [TUb). our simulations sug gest that 80% o f DLA s are within 6 p h ys = 15 hjQ kpc, which 
is consistent with the results of IGardner et al.1 (2001). However if we include the fainter 



galaxies that were not resolved in the simulations by IGardner et al.1 (120011 ). then the overall 
impact parameter distribution becomes narrow, with more than 80% of DLA sight-lines 
having b phys < 5 kpc. 

The small impact parameters suggest that DLAs in our simulations are compact, and 
this has significant implications for the nature of DLAs. The compactness of the simu- 
lated DLAs can also be observed in the projected distribution of DLAs in Figure 2 & 3 of 
Nagamine et al.l (l2004d ). There, one can see that the DLAs are at the centers of galaxies 
and coincide with star-forming regions quite well, covering roughly half the area of the star- 
forming region. These results are at odds with the idea that DLAs mainly originate from 
gas in large galactic disks (e.g., IWolfe et al.lll986l ; iTurnshek et al.lll989l ; iProchaska fc Wolfe 
19971. ll998T). T h e overa ll agreement between the simulations and observations found by 
Nagamine et al. ( 2004bl lcT) was encouraging, but the simulations with galactic wind feed- 
back now seem to underpredic t the DLA rate-o f -incid ence (see Figure [8]) compared to the 
new observational estimate by IProchaska et al.l (120051 ) utilizing the SDSS Data Release 3. 
This is because our simulations underpredict the column density distribution /(A/hi) at 
20 < log A/hi < 21 by a factor of 2 — 3 even in our h ighest resolution run with a smooth- 
ing length of comoving 1.23 h~ x kpc. ICen et al.l (120031 ) reported that they did not have this 
problem, but instead they significantly overpredicted /(A/hi) at log A/hi > 22 as well as f2 gas , 
and had to introduce dust extinction to match observations. The inconsistency in /(A/hi) 
between SPH simulations and observations and the compactne ss of the simulated D LAs 
might be related to the well-known angular momentum problem ([Robertson et al.ll2004l and 
references therein) in hydrodynamic simulations, where simulated disk galaxies are known 
to be too small and centrally concentrated. This is a problem that is not easy to solve, and 
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future studies are needed using both Eulerian mesh simulations and SPH simulations. In 
summary, there persist tensions between the predicted and observed /(A/hi), large observed 
vel ocity widths (which favors larg e and thick disks), and the predicted compactness of DLAs 
(cf. iJedamzik fc Prochaskalll998f ). 



8.4. On the Resolution Effects and Numerical Convergence 



Finally let us discuss the resolution and convergence issues of the present work. In 
Figures [6] and El for example, the results of the Q3 and Q5 run do not seem to have converged. 
This does not mean that our models for star formation and feedback are ill-defined. In fact, 



Springel fc Hernquistl (j2003al ) and ISpringel fc Hernquistl (j2003bl ) showed that the code and 
the simulations converge well in runs with star formation and wind feedback. The reason 
for the non- convergence seen in Figures [6] and [8] can be primarily ascribed to the restricted 
dynamic range of current cosmological simulations, which is present in almost all modern 
numerical work. 

In cosmological simulations, the issue of establishing numerical convergence is clouded 
by the fact that an increase of the resolution always modifies the problem one solves - one 
suddenly sees a whole new generation of low-mass galaxies that weren't there previously. 
This effect disappears only once the simulations reso lve all halos capable of c ooling, and it 
was achieved by the 'R-series' in the earlier work by [Nagamine et al.l (j2004bl ). in which we 
deduced our lower limit for the DLA halo mass. When lower mass halos are newly resolved 
in the Q5 run which were not present in the Q3 run, it steepens the slope of the relation 
between DLA cross-section and halo mass, reducing the relative contribution by the low-mass 
halos to the total DLA incidence. This change causes the difference between the results of 
Q3 and Q5 run in Figures [6] and EJ an d it can be regarded as a cautionary sign for the 
effect of resolution on these calculations. In the Q5 run, there are 150 dark matter particles 
for a halo with M halo = 10 8,5 h~ l M®, so the results of the Q5 run should be close to the 
convergence if not fully converged. 

Note however that we nevertheless find it highly useful to explore the consequences and 
predictions of our particular wind feedback model for DLAs, even if the present work eventu- 
ally can only be used to demonstrate the failure of the model. The model is numerically well 
posed, and therefore allows va ri ous predictions on the properties on DLAs and gala xies to be 
made jNagamine et all l2004bl lcl: INagamine et al.l 120051 : iNagamine et al.l l2004efl al l2005bl Jah . 



We are thus able to determine where our numerical model agrees and disagrees with the 
observational data on DLAs and galaxies. It is exactly such a comparison that makes DLA 
observations so useful for theory, because they can teach us in this way about galaxy forma- 
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tion. In the future when the computers become fast enough to achieve full dynamic range of 
mass and spatial scales without employing a series of simulations, we will be able to better 
address the issues of numerical resolution and convergence. 

Acknowledgements 

We thank Jeff Cooke and Jason Prochaska for useful discussions. This work was sup- 
ported in part by NSF grants ACI 96-19019, AST 00-71019, AST 02-06299, and AST 03- 
07690, and NASA ATP grants NAG5-12140, NAG5-13292, and NAG5-13381. The simula- 
tions were performed at the Center for Parallel Astrophysical Computing at the Harvard- 
Smithsonian Center for Astrophysics. 

REFERENCES 

Adelberger, K. L. & Steidel, C. C. 2000, ApJ, 544, 218 

Adelberger, K. L., Steidel, C. C, Shapley, A. E., & Pettini, M. 2003, ApJ, 584, 45 

Bouche, N., Gardner, J. P., Weinberg, D. H., Dave, R., & Lowenthal, J. D. 2005, ApJ, 628, 
89 

Bouche, N. & Lowenthal, J. D. 2003, ApJ, 596, 810 
-. 2004, ApJ, 609, 513 

Bouwens, R. J., Blakeslee, J. P., Illingworth, G. D., Broadhurst, T. J., & Franx, M. 2004, 
ApJ, 611, LI 

Bouwens, R. J., Broadhurst, T. J., & Illingworth, G. D. 2003, ApJ, 593, 640 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Cen, R., Nagamine, K., & Ostriker, J. P. 2005, ApJ, 635, 86 

Cen, R., Ostriker, J. P., Prochaska, J. X., & Wolfe, A. M. 2003, ApJ, 598, 741 

Chen, H.-W. & Lanzetta, K. M. 2003, ApJ, 597, 706 

Cooke, J., Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2006, ApJ, 636, L9 
Dave, R., Hernquist, L., Katz, N., & Weinberg, D. H. 1999, ApJ, 511, 521 



-22 - 



Davis, M., Efstathiou, G., Frenk, C. S., k White, S. D. M. 1985, ApJ, 292, 371 
Eisenstein, D. k Hu, P. 1999, ApJ, 511, 5 

Ferguson, H. C, Dickinson, M., Giavalisco, M., Kretchmer, C, Ravindranath, S., Idzi, R., 
Taylor, E., Conselice, C. J., et al. 2004, ApJ, 600, L107 

Fynbo, J. U., M0ller, P., k Warren, S. J. 1999, MNRAS, 305, 849 

Gardner, J., Katz, N., Hernquist, L., k Weinberg, D. H. 1997a, ApJ, 484, 31 

-. 2001, ApJ, 559, 131 

Gardner, J., Katz, N., Weinberg, D. H., k Hernquist, L. 1997b, ApJ, 486, 42 

Gawiser, E., Wolfe, A. M., Prochaska, J. X., Lanzetta, K. M., Yahata, N., k Quirrenbach, 
A. 2001, ApJ, 562, 628 

Haardt, F. k Madau, P. 1996, ApJ, 461, 20 

Haehnelt, M., Steinmetz, M., k Rauch, M. 1998, ApJ, 495, 647 

-. 2000, ApJ, 534, 594 

Hopkins, A. M., Rao, S. M., k Turnshek, D. A. 2005, ApJ, 630, 108 

Jedamzik, K. k Prochaska, J. X. 1998, MNRAS, 296, 430 

Katz, N., Weinberg, D. H., k Hernquist, L. 1996a, ApJS, 105, 19 

Katz, N., Weinberg, D. H., Hernquist, L., k Miralda-Escude, J. 1996b, ApJ, 457, L57 

Kauffmann, G. 1996, MNRAS, 281, 475 

Kulkarni, V. P., Hill, J. M., Schneider, G., Weymann, R. J., Storrie-Lombardi, L. J., Rieke, 
M., Thompson, R. L, k Jannuzi, B. T. 2000, ApJ, 536, 36 

-. 2001, ApJ, 551, 37 

Lanfranchi, G. A. k Friaca, A. C. S. 2003, MNRAS, 343, 481 

Lanzetta, K. M., Wolfe, A. M., k Turnshek, D. A. 1995, ApJ, 440, 435 

Le Brun, V., Bergeron, J., Boisse, P., k Deharveng, J. M. 1997, A&A, 321, 733 

Mailer, A., Prochaska, J. X., Somerville, R. S., k Primack, J. 2001, ApJ, 326, 1475 



-23- 

McCarthy, P. J., van Breugel, W., & Heckman, T. 1987, AJ, 93, 264 
Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
M0ller, P., Fynbo, J. U., & Fall, S. M. 2004, A&A, 422, L33 

M0ller, P., Warren, S. J., Fall, S. M., Fynbo, J. U., & Jakobsen, P. 2002, ApJ, 574, 51 
Nagamine, K., Cen, R., Hernquist, L., Ostriker, J. P., & Springel, V. 2004a, ApJ, 610, 45 
-. 2005a, ApJ, 627, 608 
-. 2005b, ApJ, 618, 23 

Nagamine, K., Springel, V., & Hernquist, L. 2004b, MNRAS, 348, 421 
-. 2004c, MNRAS, 348, 435 

Nagamine, K., Springel, V., & Hernquist, L. 2005, in IAU Symposium, ed. M. Colless, 
L. Staveley-Smith, & R. A. Stathakis, 266 

Nagamine, K., Springel, V., Hernquist, L., & Machacek, M. 2004d, MNRAS, 355, 638 

-. 2004e, MNRAS, 350, 385 

Night, C, Nagamine, K., Springel, V., & Hernquist, L. 2006, MNRAS, 366, 705 
Okoshi, K. & Nagashima, M. 2005, ApJ, 623, 99 

Okoshi, K., Nagashima, M., Gouda, N., & Yoshioka, S. 2004, ApJ, 603, 12 

Peroux, C, McMahon, R. G., Storrie-Lombardi, L. J., & Irwin, M. J. 2003, MNRAS, 346, 
1103 

Pettini, M., Shapley, A. E., Steidel, C. C, Cuby, J.-G., Dickinson, M., Moorwood, A. F. M., 
Adelberger, K. L., & Giavalisco, M. 2001, ApJ, 554, 981 

Prochaska, J. X., Herbert-Fort, S., & Wolfe, A. M. 2005, ApJ, 635, 123 

Prochaska, J. X. & Wolfe, A. M. 1997, ApJ, 487, 73 

-. 1998, ApJ, 507, 113 

-. 2002, ApJ, 566, 68 

Prochaska, J. X., Wolfe, A. M., Tytler, D., Buries, S., et al. 2001, ApJS, 137, 21 



Rao, S. M. & Turnshek, D. A. 2000, ApJS, 130, 1 

Robertson, B. E., Yoshida, N., Springel, V., & Hernquist, L. 2004, ApJ, 606, 32 
Schaye, J. 2001, ApJ, 559, LI 

Shapley, A. E., Steidel, C. C, Adelberger, K. L., Dickinson, M., Giavalisco, M., & Pettini, 
M. 2001, ApJ, 562, 95 

Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 119 

Springel, V. 2005, MNRAS, 364, 1105 

Springel, V. & Hernquist, L. 2002, MNRAS, 333, 649 

-. 2003a, MNRAS, 339, 289 

-. 2003b, MNRAS, 339, 312 

Storrie-Lombardi, L. J. & Wolfe, A. M. 2000, ApJ, 543, 552 

Theuns, T., Viel, M., Kay, S., Schaye, J., Carswell, R. F., & Tzanavaris, P. 2002, ApJ, 578, 
L5 

Turnshek, D. A., Wolfe, A. M., Lanzetta, K. M., Briggs, F. H., Cohen, R. D., Foltz, C. B., 
Smith, H. E., & Wilkes, B. J. 1989, ApJ, 344, 567 

Weatherley, S. J., Warren, S. J., Mdller, P., Fall, S. M., Fynbo, J. U., & Croom, S. M. 2005, 
MNRAS, 358, 985 

Wolfe, A. M. & Chen, H.-W. 2006, ApJ, 652, 981 

Wolfe, A. M., Lanzetta, K. M., & Foltz, C. B. 1995, ApJ, 454, 698 

Wolfe, A. M., Turnshek, D. A., Smith, H. E., & Cohen, R. D. 1986, ApJS, 61, 249 



This preprint was prepared with the AAS IATgX macros v5.2. 



-25 - 



Table 1. Simulations 



Run 






^gas 


e 


wind 


03 


144 3 


2.42 x 10 7 


3.72 x 10 6 


2.78 


none 


P3 


144 3 


2.42 x 10 7 


3.72 x 10 6 


2.78 


weak 


Q3 


144 3 


2.42 x 10 7 


3.72 x 10 6 


2.78 


strong 


Q4 


216 3 


7.16 x 10 6 


1.10 x 10 6 


1.85 


strong 


Q5 


324 3 


2.12 x 10 6 


3.26 x 10 5 


1.23 


strong 



Note. - Simulations employed in this study. All 
simulations have a comoving box size of 10 h^ 1 Mpc. 
The (initial) number of gas particles Np is equal to 
the number of dark matter particles, so the total par- 
ticle count is twice N p . m DM and m gas are the masses 
of dark matter and gas particles in units of /i _1 M Q , 
respectively, and e is the comoving gravitational soft- 
ening length in units of h" 1 kpc, The value of e is a 
measure of spatial resolution. 



Table 2. Test runs with different wind parameters 



Run 


Cdccoup 


^-maxtravel \h kpc] 




Q3 


0.1 


20 


default values 


Q3w0 


1.0 


20 


higher p decoup 


Q3wl 


0.01 


20 


lower pdccoup 


Q3w2 


0.1 


4 


shorter £™* up 


Q3w3 


0.1 


100 


longer t^ up 



Note. - - The density p d ccou P = Cdccoup Pth denotes 
the threshold density for the decoupling of the hy- 
drodynamic force, and t™^ up = Caxtravei/^ param- 
eterize the maximum allowed time of the decoupling. 
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Table 3. Fit parameters for test runs 



Run 


slope a 


normalization (3 


Q3 


0.84 


3.98 


Q3w0 


0.83 


3.96 


Q3wl 


0.80 


3.99 


Q3w2 


0.80 


3.95 


Q3w3 


0.82 


3.97 



Note. - Fit parameters to the 
DLA area vs. halo mass relation- 
ship. 



Table 4. DLA halo masses 



Run 


log M peak 


logM 50 


logM 75 


log(M DLA ) 


(logM DLA > 


03 


(8.5) 


10.1 


11.1 


11.8 


10.4 


P3 


9.6 


10.4 


11.5 


11.9 


10.6 


Q3 


10.8 


10.7 


11.7 


12.0 


10.8 


Q4 


11.6 


11.1 


12.0 


12.2 


11.1 


Q5 


12.0 


11.5 


12.3 


12.4 


11.3 



Note. - - Various DLA halo masses in units of h l M Q that 
characterize the differential distribution dN '/ '(dz d\og M) with 
the distribution cutoff at M halo = 10 85 h' l M & . The following 
quantities are listed: peak halo mass M pea k, median halo mass 
M 50 (i.e., 50 percentile of the distribution), 75 percentile of the 
distribution M 75 , mean halo mass (M DL a) (Equation (fTUj) ). and 
the mean of log Mh a ,i n rather tha n Mh a ,i n itself for the comparison 
with the result by lBouche et all (120051 ). The peak halo mass for 
03 run is shown in the parenthesis because it is the cutoff mass 
itself of the distribution. 
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Fig. 1. — Comparison of A^hi [cm -2 ] (top two panels and the range of values shown on the 
side with the color gradient) and DLA (bottom two panels) distribution in a dark matter 
halo of a mass ~ 1.7 x 10 12 h~ l M Q at z = 3 for the 03 (no wind) and Q3 (strong wind) 
run. Each postage stamp is about comoving AOOh^ 1 kpc (physical ~ 100/i _1 kpc) on a side. 
The DLA cross section in the 03 run is slightly larger than that of the Q3 run. 
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i ■ • r- 
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Fig. 2. — Distributions of Nm (top row) and DLAs (bottom row) for the same dark matter 
halo as in Fig. [T] in the four test runs with the same color scale. Each postage stamp is about 
comoving 400/i _1 kpc (physical ~ 100/i~ 1 kpc) on a side. 
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Fig. 3. — DLA cross section vs. halo mass at z = 3 for the four test runs with different wind 
parameters. Open triangles are the median in each halo mass bin. The red solid line is the 
power-law fit to these median points (Table [3]), and the black dashed line is the fit to the 
original Q3 run. The shaded contour in the background is the actual distribution of halos 
equally spaced in logarithmic scale. 
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Fig. 4. — Panel (a): Hi column density distribution function at z = 3 for the four test runs 
with different w i nd pa rameters. The data points and the red short-dashed line are from 
Prochaska et al.l (120051 ) . The solid and open triangle is the results from the original Q3 and 
Q5 runs. Panel (b): Cumulative rate-of-incidence as a function of IZab magnitude for the 
test runs and the original Q3 run. The results of the test runs are not so different from the 
original Q3 run. 
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Fig. 5. — Dark matter halo mass function of Sheth fc Tormenl ( 1999 ) at z = 3, shown in the 
form of M ■ n(M) to ease the interpretation of Equation (HI). The red dash e d line shows a 
power- law n(M) oc M~ 2 . We use the transfer function of lEisenstein &: Hul (1l999l ) and the 
power spectrum is normalized to as = 0.9. 
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log V c [km/s] 




9 10 11 12 13 14 

log M halo [h-i M ] 

Fig. 6. — Differential distribution of DLA rate-of-incidence dN/ (dz d log M) at 
function of halo mass, computed using Equations ([2D and PJ. The open crosses indicate the 
median halo masses M50, and the open triangles indicate the 75 percentile of the distribution 
M75. The filled squares indicate the mean DLA halo masses (Mdla)- The top axis also 
indicates the circular velocity as computed by Equation fJTT]) . The differences in the shape 
of the distribution arise from the differences in the relationships between DLA cross section 
and halo mass as given by Equation (J3j). 
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Fig. 7. — Halo mass vs. absolute AB magnitude at 1655A at z — 3. See text for the 
description of the three dashed lines. 
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Fig. 8. — Cumulative distribution function of DLA rate-of- incidence as a function of apparent 
IZab magnitude. The spectroscopic limit TZab = 25.5 and the limiting magnitude TZab = 30 
are indicated by the vertical dotted lines. The top axis gives the corresponding halo masses 
using the relationship for the Q3 and Q5 runs; IZab = — 2.51ogMy ia .i r , + 57.03 — 5 log h 70 



(Eq. [15]). The yellow band is the observational estimate by lProchaska et al.l (12005 
for details). 



see text 
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Fig. 9. — Cumulative probability distribution of DLA rate-of-incidence as a function of 
apparent IZab magnitude at z — 3. Panel (a) is when the distribution is normalized by the 
value at M h£do = 10 9 - 8 h^M®, i.e., assuming there wou ld be no PL As i n halo s less massive 
than this value. This roughly reproduces the result of lHaehnelt et al.l (120001 ) . Panel (b) is 
when the distribution is normalized by the value at Mh a io = 1O 8 ' 5 /i _1 M , which is more 
consistent with the results of our SPH simulations. The top axis shows the scaling with halo 
mass for the Q3 and Q5 runs, TZab = — 2.51ogM ha i + 57.03 — 5 log/i 70 (Eq. [T5].) 



- 36 - 



0.8 



X! 

A 

(X 



0.6 



0.4 



0.2 



0!5 



a 



[arcsec] 

— i — i — i — i — i — i — I — i — i — i — i— 
i l!5 

03 (no wind) 

03 (R AB <30) 
03 (R AB <28) 

P3 (weak wind) 

Q3 (strong wind) 



I 



\ \\ ■ N. 



0!5 




i i i I i i i i T~ 
10 15 

b phys i h li k P C ] 

[arcsec] 

' i ' ' ' 'l! 5 ' ' ' ' I 

— Q3 (all) 
Q 3 (R AB <30) 

Q5 (all) 

-Q5 (R AB <30) 



5 10 

b phys i h 70 k P C ] 



15 



Fig. 10. — Cumulative probability distribution of DLA rate-of-incidence as a function of 
projected impact parameter <i p h y s in units of physical hjQ kpc and angular scale in units of 
arc second at z = 3. Panel (a) shows the results of 03 (no wind), P3 (weak wind), and Q3 
(strong wind) runs. Panel (b) shows the results of Q3 and Q5 runs. In both panels, we also 
show the results when the search for the nearest galaxy is limited to those brighter than a 
certain TZab magnitude. 
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Fig. 11. — Cumulative, co moving number density of DLAs as a function of apparent TZ ab 
magnitude for 3 different values of assumed DLA physical radius r DL A = 1 kpc, 5 h^ kpc, 
and 20/ifokpc. The area of each DLA is Adla = a^dla- For each value of tdla, four 
different lines are shown: Q5 run (blue solid line), 03 (black dashed), P3 (green short-dash 
long-dash), and Q3 (red long-dashed). The thick black dot-dashed line is the cumulative 
co moving number density of LB Gs obtained by integrating the observed luminosity function 
by lAdelberger fc Steidell (120001 ) . The data point at TZab = 25.5 shows the comoving number 
density of LBGs N LBG = 4 x lO" 3 ^" 3 Mpc^ 3 . 



